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Abstract 

A simple method to calculate dispersion of the total number of 
droplets appeared in the process of nucleation caused by the smooth 
variation of external conditions has been presented. The analytical 
result for dispersion is compared with results of numerical simulations 
and the coincidence has been observed. The role of stochastic ap- 
pearance of several first droplets in formation of dispersion has been 
analyzed. 

One of the evidently characteristic features of a nucleation process is an 
occasional appearance of embryos of a new phase. The occasional manner 
of the new embryos appearance leads to generally stochastic corrections in 
a total number of droplets Ntot- In every particular process of phase tran- 
sition the number A^^^ differs from the mean value < Ntot > averaged over 
many particular attempts of nucleation in the given system. Since the pro- 
cess of condensation is essentially non linear one the mean value < Ntot > 
which is a result of averaging of particular processes differs from the value 
Ntot TAc calculated in the theory based on averaged characteristics (TAG) 
In a free molecular regime of droplets growth considered here the number 
< Ntot > approaches to the value calculated on the base of averaged charac- 
teristics Ntot TAC- It occurs when the volume V of the system is infinitely big 
which explains notation Ntot tag =< NtotiV = oo) >. One has to take this 
property into account in all theoretical descriptions of kinetics of the global 
process of phase transition. 

Earlier stochastic effects were investigated in P, jH], [Zj, but here it will 
be shown that the consideration made in [Hj, [Z] isn't satisfactory. Still 
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the structure of analysis proposed in jB], jH], [7] can be adopted in general 
features for further considerations. From the general limit theorem for the 
chain of stochastic events it is clear that the distribution P{Ntot) has to be a 
gaussian distribution with some mean value < Ntot > and dispersion 
The task is to determine these values. 

From the first point of view D equals to the value 

Do = 2< Ntot > 

of dispersion of distribution of independent events. But the nonlinearity 
of kinetics of the nucleation process leads to the deviation of D from Dq. 
Namely this deviation will be the subject of investigation in this paper. 

Here the situation of the smooth variation of external conditions will be 
studied. Namely in this situation the method proposed in 0, [S], III leads 
to the most striking error. 

The plan of investigation will be the following 

• In the first part we shall see why the previous publications jS], 05 IZ] 
can not be applied, what are the errors made in these publications. 

• In the second part a method to determine dispersion will be presented. 
This method is based on the monodisperse approximation proposed in 
jllj . To calculate dispersion one has to take into account the specific 
internal balance property in the nucleation process. As a result a very 
simple method to calculate dispersion is given. 

• In the third part the role of the several first droplets in formation of 
dispersion is analyzed. The result of this section is a simple method 
with clear physical reasons which can be applied for practically all 
significant relative deviations of from < N >. 

1 Two cycle models of nucleation process 

The formulation of the problem is described in [HI, 0; IZ| in details and one 
can omit it here. To investigate stochastic effects in kinetics of the nucleation 
process one has to write analytical approximations to the size spectrum and 
then to variate parameters of these approximations. 

The most traditional way to get analytical approximations is the iteration 
procedure. Originally started by Kuni [9_^ it was seriously modified in j2] , j3] 
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and became a powerful tool to describe nucleation kinetics. But here one can 
not use iteration procedure due to the following reasons: 



• The simplest effect which can cause the difference D from Dq is the 
influence of stochastic deviation of the rate of formation from the mean 
value on the rate of formation of droplets in next moments of time. 

• The first iteration in the iteration procedure corresponds to the ideal 
rate of growth. So, there is no influence of the deviation of the rate of 
formation of droplets of big sizes on formation of new droplets. 

• The second iteration can not be calculated analytically. 

• The standard approximations [Sj which allow to calculate further iter- 
ations correspond to some model approximation to the size spectrum. 
In these approximations one can not effectively calculate the influence 
of stochastic deviation of the rate of formation from the mean value on 
the rate of formation of droplets in next moments of time. 

The mentioned difficulties is the reason to calculate stochastic effects on 
the base of some simple models of evolution. 

To present the models for the size spectrum one has choose a renormal- 
ization in the corresponding theory based on averaged characteristics: 

• The evolution equation for the renormalized number of molecules g in 
a new phase looks like 



is the size spectrum. 

• The characteristic time necessary for decompositions (see [2]) is 
chosen to get at 2; = the maximum if the spectrum (or the maximum 
of spectrum in the first iteration). It leads to 




Here 



(f) = exp(z — g{x)) 



A =1/6^ 0.189 
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The manifestation of stochastic effects requires to have at least two effects 
- the effect of stochastic deviations and the effect of further reaction on these 
deviations. So, the simplest model is the two cycle model where during 
the first cycle the independent stochastic deviations appear and during the 
second cycle the reaction on these deviation takes place. The "minimal" 
model is the two cycle model and now we shall investigate this model. 

At first we shall consider approximation proposed in [7] . The approxima- 
tion is the following one 

• Before z = Xb = the rate of nucleation is supposed to be the ideal 
one. 

• After z = Xb = the rate of nucleation I (in renormalized units) is 
given by 

/ = exp(2; - g-) 

where is the renormalized number of molecules in the droplets 
formed before x = Xb 

In the theory based on the averaged characteristics 
g^ = A exTp{x){z — x)^dx 

J —oo 

which explains the sense of approximation. 

It is clear that this approximation is very rough. The boundary parameter 
Xb was chosen in a bad style. The reasons of such conclusion are the following: 

1. The spectrum has very discontinuous character. The jump at Xb is 
very big. It can be seen in fig.l 

We see that the region x > is approximately negligible in the calculation 
of total number of droplets. But namely this tail will be the main source of 
the deviation of dispersion from the standard value. 

2. For arbitrary Xb one can calculate the number < NfotiV = oo) > of 
the total number of droplets. But always we get the number which is greater 
than the precise value < NtotiV = oo) >= 1.002. The figure 2 shows the 
number of droplets as function of Xb- 

One can see that the number of droplets has minimum at Xboundary = 
—0.85. This value seems to be more realistic because the jump in the spec- 
trum (this jump must be in every model of such a type) isn't so giant. The 
spectrum is shown in figure 3. 
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Figure 1: The form of approximation with a boundary value Xboundary — 

The models with two cycles of evolution where the first cycle is an ideal 
one and the second cycle is formed under the influence of the first cycle will 
be referred to as "two cycle models". 

The models where the spectrum in the first cycle is determined explic- 
itly on the base of ideal conditions of the embryos formation will be called 
"explicit two cycle models". 

From figure 2 one can see that explicit two cycle models can not reproduce 
the right number of droplets. But these models have still more bad features 
concerning with calculation of dispersion. 

2 Numerical calculation of dispersion 

The best way to check the theoretical calculation of dispersion in the pre- 
sented model is to compare it with the results of numerical simulation. Here 
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Figure 2: The total number of droplets given by the theory based on averaged 
characteristics as a function of 

the following type of numerical simulation will be presented: The evolution 
of the system of the volume V starts from big negative z {z ^ —15 is suffi- 
cient) and at every step dz the random number between and 1 is generated 
V times. Every time the generated random number is compared with f * dz 
where / is the rate of nucleation in TAG and can be calculated as 

/ = exp(^-^(^)) 

To know g one has to know four integrals Ii,i — 0,1,2,3 and g is determined 
as 

g = Iq * z^ — 3z'^ * Ii + 3z * h — h 
The integrals /j are determined in TAG as 
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Figure 3: The form of spectrum approximation with a boundary value 

•^boundary 0.85. 



In application to the random appearance of embryos it means that at every 
attempt when the random number is greater than f * dz the integrals have 
to be increased as 

Certainly, at every step z moves to z + dz. 

In the current calculations the step is dz — 0.001, at every V around 1000 
nucleation processes were performed and then the mean value and dispersion 
were calculated. 

The interesting value is the ratio D/Dq which is plotted in figure 4. Since 
there is no evidence what total number of droplets one has to take in cal- 
culation of Dq (there are two variants Dq — 2 < Ntot > and Dq — 2 < 
Ntotiy — oo) > ) two curves are presented. Since < Ntot > is very close to 
< Ntotiy = oo) > these curves are close. The jumps of curves correspond to 
stochastic nature of calculations and one can see what characteristic error is 



caused by averaging only over 1000 attempts. 
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Figure 4: Relative square of dispersion as the function of the volume in 
precise simulation. 

One can see that the value of dispersion goes with increase of to a limit 
value around 9.5 rather quickly. Really in a volume unit there will be around 

/oo 
exp(x — g{x)) ^ A = 1/6 
-CXD 

droplets. So, at 10 -f- 15 droplets the limit value is already attained. 

Now one can turn to the analytical calculations of dispersion. At first 
we shall present calculation according to formulas presented in j^. Since 
the approach in [Zj is absolutely analogous to approach jS] we don't mention 
all mistakes made in [7], [S]. They are mentioned in ^2), ^2] where these 
mistakes were corrected. 

The result of calculation is shown in figure 5. Instead of [7j, jS] we consider 
the model with an arbitrary Xb instead of = 0. Here one can see the relative 
square of dispersion as function of x^. 
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We plot the parameter e — D / Dq m. the Gaussian distribution 
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Figure 5: Relative square of dispersion as function of 

We see that the dependence of dispersion on X}, is rather essential and 
there is no evidence of the choice of x^. 

It is seen that the result of the model with x;, = is two times smaller 
than a real result. But the "optimal choice" Xf, = —0.85 gives the value 
of dispersion which is rather close to the real result. It seems that here 
the coincidence is attained. But further considerations will show that this 
coincidence in no more than occasional. 



9 



3 Disadvantages of two cycle model 



Having seen the results of two cycle models one has to mention that the cal- 
culation in these models is connected with technical difficulties and requires 
serious simplifications to come to the final result. Beside many unnecessary 
restrictions (see ^2], the linearization of the action of stochastic fiuctu- 
ations of already appeared droplets on the rate of formation of new droplets 
has been done. This linearization is necessary to perform analytical calcu- 
lations but it kills all possible deviations of the mean number of droplets 
from the value calculated in the theory based on averaged characteristics. 
So, the conclusion about the absence of this shift made in [Hj, [7] is out of 
justification, it is direct consequence of linearization. A special question why 
there is no shift in the mean number of droplets will be the matter of further 
analysis and now the dispersion with all non-linear terms of the two cycle 
model will be calculated. 

The most simple and effective way to see the non-linear corrections in 
the two cycle model is to fulfill the numerical simulation corresponding to 
this model. Then one has to change the calculation of /. Instead of / = 
exp(a; — g{x)) one has to use 

/ = exp(x) 

for X < Xb and 

/ = exp(a; - g-{x)) 

for X < Xb. The value of is calculated analogously to g but the values 
of li have to be taken here at Xb (no further increase is allowed). So, from 
the first point of view there are absolutely no difficulties in simulation of this 
model. 

The numerical simulation of explicit two cycle model shows new weak 
features of this model. 

The first principal difficulty is the divergence of the model. Really, there 
is a fixed small probability that until x = x;, no droplets at all will be formed. 
This probability can be taken from Gaussian distribution as 

P(0) ~ \ exp(- < iV,„, > /2) 

but here there is no reasons to take the Gaussian distribution valid for many 
events (here there are no events). It is more reasonable to take Poisson 
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distribution but even with this distribution one has to come to 

P(0) ~ exp(-/) 

Here / is the number of possible events until Xf,. For / one can get 

I ~ / exp{x)dx = exp(xb) 

J — oo 

So, with a fixed probability (it is very small) there will be no droplets until 
Xfe. Then g_ = and the size spectrum / = exp(x — 0) grows infinitely. 

As a result we see that in explicit two cycle models there appears a 
singularity which is certainly a main peak in integration over all possible 
situations. So, it completely destroys the possibility to get reasonable results 
in the two cycle models. 

To overcome this difficulty one has to see what happens when the number 
of droplets N^jf appeared before Xh deviates from the mean value < Nf,ff >. 
Since the problem is formulated on the level of the number of droplets^ it will 
be very profitable to express evolution in terms of the number of droplets. 
Fortunately it is possible to do. In |^ a monodisperse approximation for 
the evolution in the theory based on the averaged characteristics has been 
presented for the situation of the metastable phase decay. In [TT] the same 
task was solved for the smooth variation of external conditions. Now we shall 
use results of [llj for our purposes. 

Briefiy speaking one can say that according to results of [TTj during the 
nucleation period the vapor is consumed in the unit volume by A^o =< 
Neff{V = oo) >= 1/27 droplets and this value in renormalized values is 
one and the same. 

What will happened when aX z = —3 (this is the average value for the 
position of formation of the maximum of droplets) the number of appeared 
droplets will be less than 1/27? The system will wait until the moment when 
there will be 1/27 droplets in the unit of volume. It will occur at some 
moment later than z = —3, i.e. at z = —3 + 6. Later the evolution will be 
the same but in the shifted z + 6. 

We see that if we calculate the number of droplets precisely at x = Xb then 
the parameter in a functional approximation for the rest of size spectrum 

^It is really the most important characteristic, but generally speaking one can investi- 
gate some other characteristic, for example /i(oo) — J__^ (ia;exp(a; — g{x)). 
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will be smaller or greater than the necessary value. In the monodisperse 
approximation this functional approximation is 

/ = exp(x - No{x - XbY) 

In explicit two cycle boundary it is 

f = exp{x - g_{x)) (1) 

When we put the smaller A^o than 1/27 then the validity of approximation 
will be violated. One simply can not use approximation when the value of 
parameter seriously differs from the value 1/27. 

The same is valid also for approximation (^. Approximation can 
be used only if all Jj are near J^''^x'' exp{x)dx. The use of (|T]) for all devi- 
ations (also corresponding to the small values if Jj) is the serious mistake. 
Namely this mistake doesn't allow to use the two cycle models with a fixed 
boundaries. 

To formulate the two cycle models with floating boundary one has to 
formulate the balancing property of the kinetics of nucleation. 

4 Balance property in kinetics of nucleation 

The calculations of stochastic effects are so complex that they can be made 
only in some primitive models of evolution in frames of mean values of char- 
acteristics. So, we have to present simple and rather precise models of evo- 
lution. 

At first we consider the situation of decay. In the situation of decay 
there are two sources of balancing force. The first one is evident. Suppose 
that imaginary the whole period is divided into two arbitrary parts. Let the 
number of droplets in the first part be greater that the average value. Then 
these droplets will consume metastable phase with greater intensity and the 
second part of nucleation will be shorter. Then the supposed number of 
droplets born in the second part will be smaller. This balances the number 
of droplets near the average value. 

The second balancing reason is more specific. Suppose that the first 
droplet isn't going to appear. What will be changed in kinetics? Noth- 
ing. The start of nucleation process is the moment of appearance of the 
first droplet and the system will simply wait until this moment without any 
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changes. The process of nucleation is invariant to the starting point which 
is the appearance of the first droplet. 

Denote as tgt the moment of appearance of the first droplet. Then 

< dNtot > ^ Q 
dtst 

in the situation of decay. 

Consider now the situation of smooth variation of external conditions. To 
get the analytic approximations one has to choose some model of kinetics. 
The model of monodisperse spectrum is rather simple and precise We 
shall follow this model. 

Certainly the first reason of balancing effect takes place also in dynamic 
conditions. But the second reason doesn't take place. The problem is that 
there is no fixed initial point such as tgt in decay. 

Ordinary the coordinate of monodisperse spectrum is^ z = —3. We shall 
call this coordinate z^. Then the distribution function has the form 

f = f,exp{x-N,ffix + 3f) 

where the effective number of droplets N^fj is given by 

Neff = ^, c= 0.189 

or 

Neff = — 

en 27 

The value Neff has the gaussian distribution 

PiNeff) ~ exp( 2 < Neff > ^ 

where < N^ff > is the mean value. 

Consider what will happen when aX z = —3 the required number of 
droplets N^ff will be less than^ < N^ff >. The system will simply wait until 
< Neff > will be attained. And then the evolution will be absolutely the 

^In free molecular regime when the distribution function has a maximum at z = 0. 
■^Certainly, it is simplification to grasp the main features. 
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same. But the evolution will occur with a shifted value of argument z which 
corresponds to the shift in the amplitude value 

/* exp{kSz)f^ 

where 6z is the shift in x and k is the coefficient, which can be put to 1 by 
appropriate scale of variables z, x. We keep here k to show that the effects 
don't occur due to the special renormalization. 

Then the mean total number of droplets < Ntot > will be shifted as 

< Ntot >^< Ntot > exp{kSz) 

Our task is to obtain the shift Sz or the distribution of shifts Il{Sz). This 
distribution is connected with P{Neff). 
One can present < iVg// > as 

exp{kx — g{x))dx 

-oo 

or approximately 



Then 



exp{kx)dx 

-oo 

dNeff 



Then 



dzh 



dN if 

P{N,ff) = n{5z)C-^)-^ = {UeMkz,)r^n{5z) 

One can also write the last relation as 

^ ^ e// 2 < N^ff > ' 

Then one can speak about the density distribution "E^Ntot) of the total 
number of droplets because Ntot =< Ntot > exp{k6z) Then 

E{Ntot)^Ii{6z)keMkSz) 



f^exp{kzb) (2) 



or 
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This leads to 

and there is no effect of shift. The distribution of the total number of droplets 
is a gaussian distribution with the ordinary mean value. 

The last relation leads to to the absence of shift effect also in the case of 
dynamic conditions. 

Certainly, one can mention that the absence of shift takes place only 
approximately. Really, having used in the same parameter k we supposed 
that the difference between ideal supersaturation and real supersaturation is 
absent. But the difference isn't absent, it is only small. Then instead of (j2I) 
we come to 

The sense of parameters is explained in For us now it is only important 
that at z = Zb = —3 the coefficient k will be approximately 

k' = k{l - exp(-3)) 

Then 

E{Ntot) ~ PiN.ff) expikzbil - (1 - exp(-3)))) 

or 

Because exp(— 3) is very small one can not see this effect clearly in ex- 
periment or in comparison with numerical simulation. 

In fact the existence of the last effect is rather doubtful and it is more 
preferable to speak about the absence of of the regular shift in the considered 
situation. 

One has to note that the property of the absence of the regular shift is very 
important in the context of stability of kinetics. Really all proposed models 
have to demonstrate the property of stability in respect to the stochastic 
perturbations. Namely the absence of the regular shift ensures this stability. 

One can also see that the property of regular shift will take place in all 
two-stage models. The two stage models of evolution are the models where 
the first part of spectrum is formed under the ideal supersaturation (it can 
be the monodisperse peak or explicitly the spectrum like exp(a;)) and the 
rest part of spectrum is formed only under the external influence and the 
influence of droplets from the first part of the size spectrum. In details this 
property will be analyzed separately. 
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5 Two cycle models with floating boundary 



The direct result of the previous balance property is the absence of correction 
terms to the mean number of droplets formed under the dynamic conditions. 
So, the problem to determine the mean number of droplets is solved. But 
dispersion remains still undetermined. 

But now one can formulate the true two cycle models. They must be 
the models with a floating boundary. The floating boundary is defined as 
a coordinate or the time moment Zf when the parameters of approximation 
concerning the first cycle attain their values used in approximation in TAG. 
In the monodisperse approximation there will be no difficulties - parameter 
Xf in every particular attempt is the coordinate when the number of appeared 
droplets attains 1/27. 

Now we shall perform simulations for explicit two cycle models with 
boundaries Xh = and Xb = —0.85. Results are shown in figure 6 and 
figure 7. 

One can see that the result of the model with Xf, = is satisfactory but 
there are absolutely no reasons to choose Xb = except historical ones [7]. 
The coincidence is occasional. The most reasonable model with Xb = —0.85 
gives a wrong result. 

In explicit two cycle models there will be a problem because different 
= 0,1,2,3 can attain necessary values exp(a;), = 0; — 0.85 in 

different moments. This is an evident disadvantage of explicit two cycle 
models. Having noticed the weak features of explicit two cycle models men- 
tioned at the beginning, we see that one has to consider the monodisperse 
model with the floating boundary. 

The model is the following: 

• The effective number of droplets is A^^q = 6/27 (the number 6 appears 
from A = 6 in the evolution equation in TAG). 

• Since the spectrum is formed in an internal point of a nucleation period 
one has some reasons to suppose that the subintegral function for g in 
a current moment of time is rather symmetrical being related to the 
maximum. Then the number of droplets which has to be formed before 
maximum has to be taken as No/2. 

• The moment Zf of the monodisperse peak appearance is characterized 
by formation of A"o/2 droplets. 
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Figure 6: Relative square of dispersion as a function of volume in the explicit 
two cycle model with a floating boundary — 0. 

• After Zf the spectrum is 

/ = exp(x - No{x - Xff) 

The numerical simulation requires to specify the procedure of calculation: 

• Until /o > No/2 = 3/27 the function / is exp(x) and stochastic ap- 
pearance of droplets is initiated. 

• When for the first time Iq > No/2 the coordinate z is referred^ as Zf. 

• In further evolution 

/ = exp(a; — Nf{x — Xf)^) 

^Whcn instead of the mentioned choice of z/ we choose z/ as coordinate when Iq attains 
No then the result of numerical simulation of this model differs from result of numerical 
simulation of the precise model. 
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Figure 7: Relative square of dispersion as a function of volume in the explicit 
two cycle model with a floating boundary — —0.85. 

where 

Nf^2*Io 

Dispersion in this model is shown in figure 8. 

One can see that the coincidence between precise simulation and simula- 
tion in monodisperse model is satisfactory. 

This approximation contains the precise number of effective droplets 
which doesn't precisely coincide with 1/27. So, here the number of effective 
droplets isn't the absolutely constant value prescribed by analytical recipes. 
But it is better to follow this procedure than to threw out the fractional 
part of the number of droplets at xj. Then even the number of droplets has 
the characteristic jumps at the values of volume n/27, where n is a natural 
number. These jumps a rather essential at moderate volumes. It is evident 
that dispersion near these values has to be greater than the ordinary value. 
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Figure 8: Relative square of dispersion as the function of volume of the 
system in the monodisperse approach. 



6 Analytical caluclation of dispersion 

The last model allows very simple calculations to get the dispersion of the 
total number of droplets. Really, the dispersion of the first part of droplets, 
i.e. the droplets which are in the monodisperse peak is zero 



since the floating boundary has to be used. The monodisperse spectrum 
absolutely governs further evolution. It is very easy to calculate the disper- 
sion of external influence on further evolution,. Really according to balance 
property and the self similarity of conditions of nucleation during the nucle- 
ation period (see ^T]) the relative deviation of the number of rest droplets 
equals to the ideal dispersion of the distribution of independent < N^fj > /2 
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droplets, i.e. 



Dexi = 2[< Nrest > /{N,ff/2)]* < N^est > (3) 

where N^est is the number of remaining droplets 

< Nrest > = < Ntot > - < Neff > 

Certainly one has to take into account that under the fluctuating conditions 
the droplets are formed stochastically with internal fluctuation characteristic 
to the independent formation of (< Ntot > — < Neff >) droplets 

Dint = 2 < Nrest > 

Then 

Dtot = Dext + Dint 

Having put < N^ff >= 6/27, < Ntot >= 1 one can come to 

Dtot ~ e2 < Ntot > 
e = 7 

The last value isn't too far from results of simulation. 

If one shall take into account the similarity of conditions of nucleation 
fn\ then it will be reasonable to put < Nrest >= 1 (the total number of 
droplets) and then 

e = 10 

which is practically the same result as in numerical simulations. 

The advantage of this method is that it can be used as the base for 
further more detailed constructions taking into account some fine features of 
nucleation kinetics. 

7 The role of first droplets in formation of 
dispersion 

The effect of the influence of several first droplets on the deviation of the mean 
value of the droplets number from the value prescribed by the theory based 
on the averaged characteristics manifests itself in the fact that corrections 
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to the mean value of the number of droplets can be completely explained by 
formation of several first (actually the very first one) droplets. Here such a 
radical result does not take place, but for a wide range of volumes of systems 
one can effectively use the explicit role of several first droplets. 

Since several first droplets (at least those in the monodisperse peak) are 
formed independently, the partial distribution of n droplets is the Poisson's 
distribution 

Pn ~ rexp(-/)/n! 
Here I is the number of possible events. Then at the initial period 

/ = exp(a;) 

and 

Pn{l)dl = Pn{x)dx 

It leads to 

Pn{x) = exp(nx) exp(— exp(x)) exp(x) 

Certainly, the main influence appears due to the first droplet. The partial 
distribution for the waiting of the first droplet is 

Pq = exp(x — exp(x)) 

and it is similar to the universal distribution observed in [3]. 

The values of dispersions at moderate V initiated by the first droplet, 
by two first droplets and by three first droplets are shown in the figure 9. 
The increase of the number of droplets taken into account corresponds to 
the increase of dispersion - the lowest curve corresponds to account of only 
the first droplet, the highest curve corresponds to three droplets taken into 
account. Here the value of e referred to the number of droplets prescribed 
by the theory based on averaged characteristics is drawn. 

It is seen that the increase in e initiated by every new curve is rather small 
and one can state that the main influence is initiated by the first droplet. 

One can see that the account only of three first droplets ensures the 
correct description of the monodisperse peak for situations where Neff < 3. 
It corresponds to the volume V3 = 3 * 27 = 81 which is already sufficient for 
majority of real practical situation where the stochastic effects can be really 
observed in practice. 
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Figure 9: Relative square of dispersion as the function of volume initiated 
by several first droplets. 

8 Concluding remarks 

The method proposed above is based on monodisperse approximation. This 
approximation is necessary because it allows to close the system of evolution 
equations on a level of the number of droplets. The use of explicit models 
with more complicated forms of spectrum used in the expression for the 
quantity of new liquid phase in corresponding TAG requires to write more 
complicated equations for other stochastic characteristics of these spectrums, 
which makes the problem rather more complex. 

Meanwhile the applicability of the monodispccc approximation is based 
on the sharpness of subintegral function in expression for the quantity of 
new phase. In the free molecular regime of growth this function looks like 
{z — x)^exp(x) in TAG and allows the monodisperse approximation. In 
the diffusion regime of growth the subintegral function in TAG will be hke 
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(z — exp(a;) and can not be approximated by a monodeisperce peak so 
effectively as in the case of a free molecular regime of growth. So, there ap- 
peared some difficulties but even in this situation the predictions of monodis- 
perse model qualitatively coincides with results of numerical simulation. 

In the diffusion regime of growth one has to put the position of monodis- 
perse peak at z — —3/2 and then the value of necessary droplets in the 
monodisperse peak will be N^/f = 2 which is approximately the same as the 
total precise number of droplets in TAG Ntac = 1-96. Then the total number 
of droplets appeared in monodisperse approximation will be Ntot mono = 1-98. 
Here Ntot mono was calculated by the following way: at z — —3/2 the num- 
ber of droplets will be the half of Nf.ff and later the number of droplets is 
calculated as 

mono 

J 3/2 

As the result we see that 

Ntac ~ A^e// ~ Ntot mono 

It means that practically all droplets have to be in monodisperse approxima- 
tion. The consequence of this is the absence of dispersion. But the property 
that all droplets are in monodisperse peak means that all of them at the 
same time govern their evolution. So, then one can attribute the dispersion 
of independent events. Then 

D^Do 

Results of simulation give 

D = (1.2 ^ 1.25)L'o 

which corresponds to above predictions. 

Now we shall refine the model. Now at 2; = —3/2 the monodisperse peak 
is formed and we calculate the number of droplets explicitly 

-3/2 



Ninit = / exp{x)dx = exp(-3/2) 
J —00 



Then the rest quantity of droplets is calculated by the previous formula. 

The result of solution is that approximately the three quarters of droplets 
have to be included in monodisperse peak and the rest quarter is formed 
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under the influence of this peak. The following question appears: shall one 
takes into account the dispersion of droplets in the peak? The previous result 
which shows the self action of the monodisperse peak onto itself requires to 
take into account the dispersion of the monodisperse peak. So we have 

D = 3/4 * + 1/4(1 + 4/3) * Do 

where Di is dispersion of initial peak. Certainly, for Di one has to take the 
dispersion of independent events Dq. In the last term 1/4(1 + 4/3) * Dq the 
value 1 comes from the self dispersion of free events and term 4/3 is the 
induced dispersion. 

Earlier in investigation of free-molecular regime there was no use to take 
into account the dispersion of initial peak because it would give addition only 
6/27. This addition is so small in comparison with renormalizing factor 10 
appeared earlier that there is no need to take it into account. In diffusion 
regime such addition can be essential. 

Then the calculation according to the last formula gives 

D = 4/3* Do 

Then the result is slightly greater than the result of simulation but the dif- 
ference is negligible. 

This result is the estimate from above because the substitution Di = Dq 
increases the result. The flrst estimate D = Dq is the estimate from below. 
So, 

Do<D < 1.33^0 

The interval ensured by these estimates is rather narrow and the result of 
simulation lies in this interval. 

One has to mention that in the diffusion regime foundations of the monodis- 
perse model and of all further constructions are rather weak. So, one can not 
state that our result isn't well justifled. Fortunately, in the diffusion regime 
of growth kinetics of nucleation is analyzed separately P and has absolutely 
another physical reasons of evolution. Also one has to take into account the 
effect of compensation due to the growing volumes considered in [T2]. 
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